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TECHNICAL FUEL© 

This invention relates to a method and apparatus that can measure the surface profile and 
twist of one or both sides of an object that is moving longitudinally relative to the apparatus, 
independent of rigid-body motions in the measurement direction. Such measurements are useful 
for surface quality control of manufactured products, for example sawn lumber and rolled metals, 
and for extended surfaces such as railroad rails and highway pavements. 

BACKGROUN© OF THE INVENTION 

Surface profiling is an important need in many industrial and scientific applications. 
Typical examples include flatness inspection of railway rails and road surfaces, and quality control 
of manufactured products such as sawn lumber and rolled metal. 

The simplest way to measure surface height profile is to make a series of measurements 
with a distance sensor while relatively moving the measured object and the sensor in a straight line 
perpendicular to the measurement direction. FIG. 1(a) schematically shows an example 
arrangement. The drawback to this method is that deviations from straight-line motion cause 
relative displacements in the measurement direction that are indistinguishable from measured 
surface shape. Thus, very accurate linear motion is essential. 

In many cases, accurate linear motion is not practicable. For example, when measuring 
the surface profile of a long length of railway track, it is not feasible to provide a separate linear 
slide for moving the sensor. Instead, the sensor must travel on the irregular track that it is 



measuring. Existing non-inertial techniques focus on the surface curvature because it can be 
identified independently of relative displacements or rotations. A typical arrangement uses three 
or more sensors operating simultaneously, as shown in FIG. 1 (b). US Patent 4,288,855 describes 
examples of this concept. The sensors estimate the local surface curvature from the second finite 
difference of their measurements. The curvature values are integrated twice to determine the 
surface profile. This method is effective, but it has difficulty resolving surface features that are 
either greatly shorter or greatly longer than the total spacing of the sensors. 

US Patent 5,280,719 describes an apparatus that uses a large number of equally spaced 
sensors. The apparatus seeks to identify long surface features by overlapping sets of 
measurements that are made at successive intervals much less than the total spacing of the 
sensors. However, the large number of sensors that are required creates a large cost and 
maintenance burden. US Patent 4,573,131 describes a method of using just four sensors to 
achieve the same objective. The sensors make successive measurements at small intervals of 
travel. The sensor spacings are chosen so that the numbers of intervals of travel between 
successive sensors have no common factors other than unity. This feature avoids the existence of 
"nulls", which is a measurement artifact whereby surface features with specific wavelengths are 
not observed. The method described in US Patent 4,573,131 is based on the assumption that 
successive sensors reaching the same point along the measured line have the same distance from 
the measured surface. This is typically not a reasonable assumption, and the described method 
has only limited effectiveness. 

An alternative approach, the "inertial" method, uses an accelerometer that runs at a 
constant speed along the surface to be profiled, giving a signal that is proportional to surface 
curvature. This signal is integrated twice to obtain the surface shape. Conceptual simplicity 
makes this an attractive method. However, there are also some practical limitations. The first 
concerns the high frequencies that need to be measured. These can extend beyond the capabilities 
of the accelerometer. A second limitation is the assumption that the accelerometer and the 
surface are permanently in contact. This is difficult to achieve reliably, especially at high scanning 
speeds due to the large inertial forces that act on the system. Finally, only smooth profiles can be 
measured. A sharp step, for example, would not be detected properly. 

None of the above techniques is well suited to measuring the surfaces of two-sided 
objects. At best, they can measure each of the two sides separately, but they do not provide 
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accurate thickness information. Similarly, they do not provide accurate twist information with 
parallel measurements, nor are they suitable for two-dimensional surface evaluations. . 

BRIEF SUMMARY OF THE INVENTION 

The present invention provides an improved method for measuring surface height profile 
that overcomes the limitations of previous art. The invention involves using multiple distance 
sensors operating simultaneously. The main difference from previous designs is the way in which 
the sensor measurements are combined to determine the desired surface height profile, and in the 
flexibility of the possible geometrical arrangements that can be used to achieve particular 
measurement objectives. 

The central idea of the present invention is the observation that surface height features 
appear in delayed sequence as the observed surface moves longitudinally relative to the sensor 
array. However, any relative motions in the measurement direction, either translational or 
rotational, appear simultaneously at all sensors. Modern mathematical procedures are used to 
separate the delayed and simultaneous components of the measurements, from which the surface 
height profile is determined. 

One of the advantages of the proposed method is that it is able to resolve both fine details 
and long wavelengths. It also enhances measurement accuracy by exploiting the geometrical 
relationships that exist among successive measurements, and by averaging the redundant 
information from multiple sensors. In addition, the proposed method is able to reconstruct the 
entire length of the scanned surface. This is in contrast to the curvature method, where the 
extremes are lost because neighboring points are used to make each curvature calculation. The 
underlying idea of the present invention is very flexible, and it can easily be adapted to 
simultaneously measuring the surfaces of two-sided objects. Further possibilities include 
measurement of surface twist and two-dimensional surface scanning. 

BRIEF DESCRIPTION OF THE DRAWINGS 
FIGURE 1. Surface profile measurement, (a) using a single sensor, (b) using multiple sensors. 
FIGURE 2. Four-sensor arrangement for one-sided measurements. 
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FIGURE 3. Distance profiles measured by sensors d-c-b-a and the calculated surface profile. 

FIGURE 4. Six-sensor arrangement for two-sided measurements. 

FIGURE 5. Six-sensor arrangement for parallel one-sided measurements. 

FIGURE 6. Eight-sensor arrangement for parallel two-sided measurements. 

FIGURE 7. Two line-type disiance sensors used to make two-dimensional surface topographic 
measurements. 



BETA1DLE© DESCRIPTION OF THE INVENTION 

To simplify the description, it is supposed that the distance sensors remain stationary and 
that the measured object moves. However, it is apparent that the description also applies to the 
case of a stationary object and moving sensors. 

FIG. 2 schematically shows an example sensor arrangement. In this case, there are four 
distance sensors arranged in a row. However, the invention is not limited to this specific 
arrangement. Other arrangements with fewer or more sensors in one or more rows, on the same 
or opposite sides of the measured object, are also useful. The distance sensors can be of any 
suitable type, chosen according to the needs of the particular application. Example sensor types 
include, but are not limited to, laser triangulation, ultrasound, eddy current, LVDT sensors, etc. 
For notational convenience, the sensors are labeled a, b, c and d, in the reverse sequential order 
that they are reached by the moving object. 

FIG. 2 schematically shows a simple roller conveyor system. This is included in the 
diagram to assist conceptual understanding, but it is not part of the present invention. Almost any 
means of longitudinally moving the measured object could be suitable. The arrangement of 
eccentrically mounted rollers in FIG. 2 symbolically indicates that relative motions in the 
measurement direction do not need to be avoided. 

The conveyor system in FIG. 2 includes a means for measuring the longitudinal motion of 
the moving object. Any accurate measuring means is acceptable; for example, a rotary encoder 
attached to a conveyor roll. The encoder readings are used to trigger simultaneous measurements 
from the distance sensors, at regular small intervals along the measured object. It is preferable 
that the sensor spacing should be an integer number of these intervals. To avoid the presence of 



"nulls" (an inability to detect surface features with specific wavelengths), the numbers of 
measurement intervals within at least two of the sensor spacings should have no common factors 
other than unity. These conditions are desirable, but not essential. Satisfactory operation can 
often be achieved even if they are not obeyed. 

As an example, consider the case of an object of total length equal to n-l measurement 
intervals (= n measured points). Let w, be the heights of the measured points on the surface of the 
object relative to a datum line fixed on or adjacent to the object. The index / is in the range 1 < * 
< w, starting at the end of the object which reaches sensor a first. Let m>, be the corresponding 
position of the datum line relative to the distance sensors. Further, let z, be the rotation of the 
datum line at the moment when point i is being measured by sensor a. Any specific center of 
rotation is potentially acceptable for defining m>, and z,. A desirable choice is the midpoint of the 
length spanned by the distance sensors because it improves the stability of the subsequent 
mathematical calculations. For dimensional consistency, the quantity z t is defined as the rotation 
angle in radians multiplied by half the length spanned by the set of distance sensors. Without any 
loss in generality, the initial quantities wj and z; can be arbitrarily assigned zero values. This 
assignment provides the definition of the position and orientation of the datum line. 

The geometry in FIG. 2 creates an inversion of the measured surface data. An increase in 
surface height causes a decrease in distance sensor reading. To avoid possible confusion, it is 
convenient to define "height" quantities of the kind: 

a, = h - a*j (1) 

where a*/ is the measured distance of the surface from sensor "a", h is a "reference" distance. 
The exact value chosen for h is not critical. For example, the average value of a*/ would be a 
reasonable choice. However, the value(s) of h chosen for the other sensor(s) cannot be chosen 
with similar freedom. They should be such that all measured "heights" would be equal for the 
case of a flat object moving parallel to the sensor array. 

For the arrangement shown in FIG. 2, the sensor measurements are related to the surface 
profile and the rigid-body motions in the measurement direction through the equations: 
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Uf + Wj + (0 - m) Zj 
U i+ P + w i + (p-m)z i /2 
"/+ + w, + (p + ? -m) z, /2 
u j+p+q+r + Wj+ (p + q + r-m)z i /2 



(2) 



where /?, # and r are the numbers of measurement intervals between distance sensors a, b, c and 
d, respectively, and m = (p+q+rJ/2. 

The relationships between the quantities in equations (2) are purely geometrical, and can 
be expressed in matrix form. For example, when n = 9 and p = 1, q = 2, r = 1: 
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corresponding to: 



(3) 



A f = g 



(4) 



where f is a vector comprising the surface profile and relative motion results, and g is a vector 
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comprising the measured distance data. Matrix A clearly demonstrates the nature of the 
measurements. The leftmost block of columns illustrates the delayed sequence of the surface 
profile data, while the center and rightmost columns of blocks show the simultaneous character of 
the relative motion data. For mathematical uniqueness, equation (3) excludes the redundant 
columns corresponding to Wj = 0 and zj = 0. 

At every measurement point, there are three unknown quantities to be determined, u iy w t 
and it. Thus, while it is possible to start measurements when the object reaches the first distance 
sensor, it is desirable to wait until at the object reaches at least three sensors. In practice, it is 
convenient to wait until the object reaches all available sensors. In this case, each of the sensors 
makes the same number of readings on the object. Under these conditions, the number of 
unknown quantities f in equation (4) equals 3n~2(p+q+r+l), and the number of measured 
quantities in g is 4(n-p-q-r). For n > 2(p+q+r-l), there is an excess of data over unknowns. This 
is desirable because the excess data can be used to reduce the effect of measurement errors. 
Under these conditions, equation (4) can be solved by the least-squares method using the "Normal 
Equations": 

A T A * f = A T g (5) 



For the example case in equation (3), the Normal Equations are: 
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For small values of n, a general matrix equation solver is sufficient to solve the Normal 
Equations. Matrix A T A is symmetrical, so some computational economy can be achieved by 
using a specialized solver for symmetric matrices. For larger values of n, iteration techniques, for 
example the Gauss-Seidel method, are useful. 

FIG. 3 shows an example measurement made using the arrangement shown in FIG. 2. In 
this example, n = 67, p = 2, q = 7, r =2. To make clear the effectiveness of the method, a small 
steps was attached to one of the rollers so that the measured object moved along an obviously 
irregular path through the roller system. The first four traces show the sensor readings and the 
fifth profile the computed surface profile. This profile corresponds accurately to the actual curved 
profile of the measured object, and successfully eliminates the effects of the rigid-body motions in 
the measurement direction. In addition, the full length of the object has been measured, even 
though each individual sensor measured less than the entire length. 

In general, it is desirable to use more than the minimum number of distance sensors that 
will give an excess of data over unknowns. The least-squares method will then have an averaging 
effect on the computed profile, and reduce the effects of random measurement errors. When high 
accuracy is required, the use of more sensors is an appropriate choice. Matrix A becomes fuller, 
but the size of A T A in equation (5) remains the same and no new ideas are required. 

When cost is a critical factor, and more modest accuracy is acceptable, it is possible to use 
fewer than the "minimum" number of sensors, for example three, or even two, instead of the four 
shown in FIG. 2. In this case, there will be fewer data than unknowns. It may also occur that for 
some reason it is not possible to avoid having common factors in at least two sensor spacings. 
Both of these cases cause equation (5) to become singular. This circumstance can be handled by 
using "regularization", a mathematical process which introduces smoothing. 

Several regularization techniques are commonly used, among them truncation and 
filtering. A convenient technique is Tikhonov regularization. The procedure involves modifying 
matrix A T A in equation (5) to penalize extreme values in the solution. The penalty can be applied 
to the norm of the model, thereby creating a "small" solution, or to the norm of the numerical first 
or second derivatives of the model, respectively creating "flat" or "smooth" solutions. The 
simplest form of Tikhonov regularization adds an extra term to the Normal Equations, as follows: 

(A T A + p B T B) f = A T g (7) 
% 



where P is a weighting factor called the regularization parameter. Matrix B is the identity, first or 
second derivative operator, chosen according to whether a "small", "flat" or "smooth" model 
solution is desired. In general, a "smooth" model is a reasonable choice because it is the most 
responsive to measurement noise. 

The Morozov discrepancy principle can be used to identify the optimal value of the 
regularization parameter p. If the standard errors, a, in the data are all equal, independent and 
Gaussian with zero mean, the optimum parameter p makes the norm of the misfit between the 
actual measured distances and the distances corresponding to the calculated solution equal to the 
standard error of the distance measurements. In this case: 



Aff - g[ ' 2 



= N (8) 



If the standard errors in the measured distances may not be assumed to be equal, then appropriate 
weighting factors are required. In that case, the Morozov discrepancy principle requires that the 
chi-squared statistic equals N. 

Given an estimate of the standard error in the measured data, o, equations (7) and (8) can 
be solved iteratively. The result is not sensitive to modest variations in P, so accurate 
convergence is not necessary. This tolerance is attractive because it could remove the need to 
evaluate p for every profile. For stable and consistent applications, it may be sufficient to use a 
fixed value of p that is updated only periodically. 

Equation (7) applies to the simplest case where only one type of regularization is used, 
and equal weighting is applied to all elements of the model vector. For specialized applications, 
equation (7) can be expanded to use all three types simultaneously, with individual matrices to 
indicate the weighting parameters for each regularization type and element in the model vector. 

It should be noted that regularization enables a surface profile to be calculated when a 
mathematical solution many not otherwise be possible. It may also be used to provide smoothing 
for a noisy solution. However, in neither case does it create any new information. Therefore, it is 
not a substitute for data from additional distance sensors. 

The idea behind this invention can be exploited further to allow efficient two-sided 
profiling using a sensor arrangement such as is shown in FIG. 4. The attraction of this 



arrangement is that the upper and lower sensors observe the same rigid-body motions as the 
single-sided arrangement shown in FIG. 2. This commonality between of both sets of sensors 
means that only one new variable is introduced, namely the lower surface profile. Thus, a 
minimum of one additional lower sensor is required to provide the needed additional data, far 
fewer than needed for two independent profile measurements. 

In the arrangement shown in FIG. 4, one sensor has been transferred from the upper 
surface, and two further sensors have been added. Although the sixth sensor is not essential, it is 
included here to provide a symmetrical arrangement with similar expected measurement accuracy 
for both sides of the measured object. The additional sensor provides redundant data that is 
averaged with all the other data. This feature enhances the measurement accuracy for both sides. 

For the arrangement in FIG. 4, the relationships of the sensor readings to the surface 
profiles and rigid-body motions are: 
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where m = (p+q)/2 = (r+s)/2. The equations involve one new surface profile variable, v /5 but no 
additional rigid-body motions. These equations can be expressed in matrix form, analogous to 
equation (3), from which the surface profiles can be calculated. As previously, at least two sensor 
spacings should not have any common factors other than unity. The skew-symmetric arrangement 
in FIG. 4 is not essential, but is chosen to avoid any left-to-right bias. An important feature of the 
use of the present invention for two-sided profiling is that it preserves geometrical thickness 
relationships, even in the presence of possible errors in identifying rigid-body motions. Thus, the 
advantages of conventional thickness measuring techniques are retained. 

Exploitation of common components in the measurements from multiple sensors can be 
carried a further step. Practical objects are often quite wide in the direction out of the paper in 
FIG. 1. In such cases, it can be of interest to measure the surface profile along parallel lines on 
one or both opposing surfaces. The availability of parallel profiles along a surface gives piore 
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than just multiple measurements of surface topography. The difference between the profiles gives 
twist, which is a geometric property of significant interest. Parallel profiling introduces a third 
significant rigid-body motion in addition to the rigid-body displacement and rotation considered 
so far, namely, rotation around the longitudinal axis of the measured object. 

Parallel profile measurements are possible for both one- and two-surface profiling. FIG. 5 
shows an arrangement of six sensors to do one-surface profiling along two parallel paths. As 
before, a symmetrical arrangement has been chosen so that equal measurement accuracy is 
expected for each measured profile. In this case, the number of sensors used equals the minimum 
required to achieve a full-rank matrix. 

For the sensor arrangement shown in FIG. 5, the relationships of the sensor readings to 
the surface profiles and rigid-body motions are: 
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where w, and v, are the two parallel surface profiles, m = (p+q)/2 = (r+s)/2, and y t = (J) A/2 
corresponds to the rigid-body rotation around the longitudinal axis of the measured object by an 
angle ((>, with a distance A between the parallel paths. As before, these equations can be expressed 
in matrix form, analogous to equation (3), from which the surface profiles can be calculated. 

FIG. 6 shows a sensor arrangement that uses eight sensors to do parallel profiling along 
four lines adjacent to the edges of a rectangular object. This arrangement combines the 
measurement features of FIG. 4 and FIG. 5 

FIG. 7 shows another useful measurement arrangement. In this example, two line-type 
distance sensors replace the single-point-type previously discussed. These sensors measure 
distance at a sequence of many points along two lines across the surface of the measured object. 
To avoid "nulls" in the subsequent computations, the lines should not be parallel to each other. 
The attractive feature of the sensor arrangement in FIG. 7 is its ability to identify a two- 



dimensional topographic map of the surface height of the measured object, independent of rigid- 
body motions. Clearly, the amount of data to be handled from the line-type distance sensors in 
FIG. 7 is much larger than from single-point-type. However, the procedures for working with 
those data still involve the same principles of instantaneous observation of rigid-body motions, 
and delayed observation of surface features. Thus, the mathematical calculations follow the same 
patterns as previously described, with some adaptations to allow for the possible non-integral 
number of measurement intervals between corresponding pixels along the measured lines. 

Commonly, the number of measurements made, using either single-point-type or line-type 
distance sensors, can be very large. In some cases, for example, when measuring road surface 
profiles, the measured object has "infinite" length, resulting in an indefinitely large number of 
measurements. In such cases, evaluation of matrix A T A for solution of equation (5) or (7) is not 
practicable. Instead, a progressive iterative solution method can be used, where the first part of 
the matrix is solved before all the remaining parts are available. 

The first step when formulating a progressive solution is to reorder the variables from the 
arrangement shown in equation (6). The original arrangement was chosen because it reveals the 
matrix structure. However, sequential arrangements of the form J = [uj wj zj u 2 w 2 z 2 . . .] T 
and g = [aj bj cj dj a 2 b 2 c 2 d 2 . . ,] T are more practical for iterative solutions. Even for 
whole-matrix solutions, the sequential arrangements are desirable because they reduce matrix 
bandwidth. 

The following is an example progressive iterative solution method for equation (6), based 
on the Gauss-Seidel method. Successive over-relaxation could additionally be used. 

Initial estimate for the first 5 data sets: 
Wj = 0 z\ = 0 U\ = ai 

For / = 2 . . 5: z, = 0 = (a 5 -di) (i-l)/4 w, = a, - w f (1 1) 

Iteration scheme: 

For j = 6 ... measure next a, bj Cj dj 
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= [dji -w H -z H + Cj-wj -.5zj]/2 

Uj+i = + c^a-HjLj -.5z ; - 2 + bj-wj+ .5 Zj] / 3 

For i = y . . . y-£ (counting in reverse) 

z f - = f -fa - ^ - .5(bi - w /+i ; + .5(d - w 5+/ ; + fa - u 4 +i) J/ 2.5 
Wj = [ fa -wj + (bi-ui+i) + (Ci-u 3 H) + (di-u 4 +i) J/4 

When 4 < / w, = / (a, - w, + z^ + fa.; - hv/ + .5z,_ 7 j + fa_ 5 - -.5z,. 5> ) + fa - w M -z^ 7 / 4 
When 3<i <4 w, = / fa - w, + z^ + fa. 7 - w lW + .5z, w > + (c^ - w,- 5 -.iz/.^ J/3 
When Ki<3 u t = f (aj - m>, + z^ + (b^ - w,.y + ,5z, w J J/2 

When /<7 skip (12) 

The variable k controls the number of iterations used to refine the solved quantities. 
Larger values give more accurate answers, but at the the expense of computational effort. A 
value of k slightly larger than p+q+r typically gives acceptable results. 

It is apparent that the ideas that constitute this invention can be applied to meet a wide 
range of measurement needs. Such needs may be met by using various numbers and types of 
distance sensors with various geometrical arrangements. The present disclosure describes the 
basic principles of the present invention, and gives many practical details. It is evident that a 
person skilled in the art could use these principles and details to construct a surface profile 
measuring device that varies in particular details from the specific examples described herein. It is 
also evident that such a person could also employ mathematical methods different from the 
examples presented here to achieve the same or similar numerical results, for example, using 
transforms. All such variations lie within the scope of this invention. 
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